function [a, b, c, absbank] = LinearBankProfile(e, e0, ef, absBank0, absBankf)
a = 0;
b = (absBank0 - absBankf)/(e0 - ef);
c = absBankf - b*ef;
absbank = a*e^2 + b*e + c;
end